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Abstract 


The  permanent  magnet  synchronous  motor  (PMSM)  is  increasingly  playing  an  important 
role  in  electric  propulsion  systems  due  to  its  many  advantages  over  competing 
technologies.  For  successful  implementation  of  the  PMSM,  rotor  position  and  speed 
information  are  required.  A  resolver  or  encoder  attached  to  the  shaft  of  the  machine 
usually  provides  this  information.  Many  applications,  however,  cannot  tolerate  the  use  of 
the  position  sensor  due  to  space  and  weight  limitations,  reliability  concerns,  or  packaging 
issues.  Thus,  there  has  been  an  intense  interest  in  the  development  of  a  so-called 
sensorless  drive,  where  the  PMSM  stator  itself  is  used  as  the  rotor  position  sensor.  In  this 
work,  a  previously  developed  sensorless  electric  drive  technique  [1]  is  applied  to  a  high 
power  permanent  magnet  machine.  The  sensorless  drive  hardware  implementation  is 
presented  in  detail.  A  set  of  experiments  is  used  to  quantify  the  performance  of  the 
sensorless  drive  approach  under  a  wide  variety  of  operating  conditions,  and 
implementation  issues  and  results  are  discussed.  The  developed  sensorless  electric  drive 
serves  as  a  technology  demonstrator  for  future  naval  applications. 
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INTRODUCTION 


The  recent  advances  in  high  power  semiconductor  switching  devices,  magnetic  materials, 
and  energy  storage  systems  tove  generated  considerable  interest  in  electric  drive  systems 
for  vehicle  propulsion  applications.  Electric  drive  systems  for  propulsion  can  eliminate 
the  need  for  shafts  and  gearing  while  increasing  vehicle  stealth  and  flexibility.  Novel 
propulsion  configurations  and  control  strategies  are  pushing  the  need  for  high  power 
density,  high  performance  electric  machines  to  operate  in  the  harsh  environments  present 
in  many  electric  drive  applications. 

The  permanent  magnet  synchronous  motor  (PMSM)  offers  several  features  that  make  it 
attractive  for  use  in  electric  drive  systems.  The  PMSM  offers  increased  power  density, 
efficiency,  and  excellent  acoustic  performance  when  conqjared  to  induction  and  DC 
motors.  In  order  to  fully  exploit  the  characteristics  of  the  PMSM,  high-resolution  rotor 
angle  information  is  required.  This  requirement  originates  from  the  need  to  effectively 
perform  torque  control  on  the  PMSM.  Furthermore,  for  speed  control  applications,  the 
speed  feedback  is  also  required.  A  resolver  or  encoder  attached  to  the  shaft  of  the 
machine  usually  provides  both  the  position  and  speed  feedback.  Many  applications, 
however,  cannot  tolerate  the  use  of  such  a  position  sensor.  Thus,  there  has  been  an 
intense  interest  in  the  development  of  the  sensorless  drive,  where  the  PMSM  stator  itself 
is  used  as  the  rotor  position  sensor. 


1.1  Statement  of  Problem 

Along  with  the  recent  interest  in  the  PMSM  by  both  military  and  industry  came  the 
development  of  new  control  techniques  for  current  and  direct  torque  control  [2,  3], 
efficiency  maximization  [4],  and  torque  ripple  minimization  [5,  6,  7,  8].  Each  of  these 
control  techniques  require  accurate  knowledge  of  the  rotor  angle  in  order  to  be  effective. 
At  the  same  time,  new  motor  propulsion  topologies  such  as  integrated  motor/propulsors 
that  place  the  motor  directly  in  a  submerg^  environment  [9]  are  initiatii^  the  need  to 
eliminate  these  rotor  position  sensors.  The  requirement  for  high  resolution  rotor  position 
information  coupled  with  the  unreliability  of  the  rotor  position  transducers  in  an 
underwater  environment  has  created  a  need  for  a  rotor  position  sensorless  strategy 
suitable  for  such  propulsion  applications.  Thus,  the  development  of  a  sensorless  drive 
strategy  for  the  PMSM  is  an  enabling  technology  toward  the  successful  implementation 
of  the  integrated  motor/propulsor.  The  elimination  of  the  high  resolution  rotor  position 
sensor  enhances  the  value  of  using  the  PMSM  in  many  other  Navy  applications  as  well 
due  to  space  and  weight  savings  that  translate  to  additional  payload  volume  and  increased 
reliability. 


1.2  Statement  of  Purpose 

In  this  dociunent,  the  implementation  of  a  sensorless  PMSM  drive  is  presented.  The 
hardware  implementation  is  detailed,  implementation  issues  are  discussed,  and 
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experimental  results  are  used  to  verify  the  effectiveness  of  the  developed  system. 
Through  experimentation,  the  developed  sensorless  strategy  will  be  shown  to  result  in 
Naval  electric  drive  capabilities  that  folly  utilize  many  of  the  attractive  characteristics  of 
the  PMSM  without  the  requirement  for  a  high  resolution  position  sensor. 


1.3  Organization  of  Report 

This  document  focuses  on  the  implementation  of  a  sensorless  electric  drive  for  a  high- 
power  PMSM.  In  Chapter  2,  a  cursory  review  of  the  rotor  position  estimation  technique 
developed  in  [1]  is  presented.  To  conqjlement  this  original  development,  a  rotor  angle 
estimation  error  a^ysis  is  performed  to  quantify  the  effects  of  PMSM  model 
uncertainty. 

In  Chapter  3,  the  hardware  implementation  of  the  rotor  position  estimation  technique  is 
presented.  The  details  of  the  hardware  implementation  are  shown  in  detail,  and 
implementation  issues  are  discussed.  The  performance  of  the  resulting  sensorless  drive  is 
examined  through  the  extensive  real-time  experiments  in  Chapter  4.  Finally,  conclusions 
of  the  research  and  remaining  issues  are  set  forth  in  Chtqjter  5. 


2  ROTOR  POSITION  ESTIMATION 

In  this  section,  a  review  of  a  model-based  rotor  position  estimator  developed  in  [1]  is 
presented.  From  this  proposed  rotor  position  estimator  model,  rotor  angle  estimation 
error  analysis  is  performed  so  that  the  effects  of  PMSM  model  uncertainty  and 
measurement  errors  can  be  quantified. 


2.1  Sensorless  PMSM  Model 


The  flux  linkage  model  of  the  non-salient  PMSM  in  the  two-phase  stationary  reference 
frame  [10]  is: 


A' 

/ 

—  r 

~ X  n'i 

^pma 

yp. 

\ 

Xp 

(1) 


where 


(2) 


The  terms  and  ^  represent  the  flux  linkages  and  voltages  corresponding  to  the 
two  fictitious  windings  in  the  stationary  reference  flume.  The  symbols  ,X„,R,  and 
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e  represent  the  self  inductance,  PM  flux  constant,  resistance,  and  rotor  angle, 
respectively. 

To  develop  a  model  for  a  rotor  position  observer,  the  vector  defined  in 

(2)  is  considered  to  be  a  disturbance  state  that  satisfies  a  known  differential  equation. 
Combining  the  disturbance  state  vector  with  the  state  variables  defined  in  (1)  yields  the 
metastate  vector,  x: 

X  =  [xi  ^2  X3  X4 

The  input  vector  is  the  applied  phase  voltage  in  the  stationary  two-phase  reference  frame 

u  =  [Mi«2r=t’aV^f ,  (4) 

and  the  output  vector,  y,  is  chosen  to  be  the  measurable  currents: 

y=lviy2r=[4»>F- 

Thus,  the  metasystem  state  and  output  equations  are  described  by 


respectively,  where 


x  =  Aa,x  +  Bu 
y  =  Cx  , 


0 

T 

0  ■ 

'1 

o' 

-T 

0 

X 

;  B  = 

0 

1 

0 

0 

-a. 

0 

0 

0 

0 

0 

0 

r 

c 

(6) 


(7) 


(8) 


The  term  represents  the  angular  velocity  of  the  rotor  shaft  in  electrical  radians  per 
second.  The  subscript  (o  associated  with  the  matrix  is  used  to  indicate  variation 
with  the  rotor  angular  velocity. 

Under  the  assumption  that  rotor  velocity  is  approximately  constant  over  the  system 
sampling  period,  the  angular  velocity  may  be  considered  to  be  a  slowly  varying  PMSM 
parameter.  This  assumption  is  conservative  given  the  large  mechanical  time  constant 
relative  to  the  fest  electrical  time  constant  of  the  typical  PMSM  propulsion  motor.  With 
this  assumption,  the  resufting  model  is  a  linear  time  varying  system,  allowing  the  use  of 
well-developed  linear  control  methods.  It  is  also  stressed  that  the  model  given  in  (6)  and 
(7)  does  not  contain  any  mechanical  variables  such  as  torque,  friction,  and  inertia  in  the 
state  equations.  This  is  advantageous,  since  torque  and  inertia  are  often  unknown  and 
vary  during  the  operation  of  the  system. 
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2.2  Rotor  Position  Observer 


With  angular  velocity  considered  to  be  a  slowly  varying  parameter,  (6)  and  (7)  may  be 
used  to  construct  a  full  order  observer  [11].  The  observer  is  used  to  estimate  the  state 
vector  X  of  (3)  from  knowledge  of  the  input  vector  u  and  a  direct  measurement  of  the 

output  vector  y.  Note  that  from  the  estimation  of  Xj  and  ,  the  rotor  position  estimate  9 
may  be  determined  from 


/a  ^ 


^  =  tan 


-1 


(9) 


Using  standard  observer  design  techniques  for  a  linear  system  [12],  the  form  of  the 
proposed  rotor  position  observer  is 


e  =  (A„-GC). 

where  e  is  the  state  estimation  error  vector 

e  =  x-i  , 

G  is  the  observer  gain  matrix,  and  i  is  the  estimated  state  vector. 


(10) 

(11) 


The  convergence  of  the  estimated  state  variables  toward  their  actual  values  is  achieved  by 
conventional  pole-placement  techniques.  The  eigenvalues  of  the  characteristic  equation 

|AI-(A^-GC)|  =  0  (12) 

are  chosen  to  have  negative  real  components  so  that  asymptotic  state  reconstruction  is 
achieved. 

Clearly,  observability  of  the  system  is  a  requirement  for  the  estimation  of  PMSM  states. 
The  condition  for  observability  is  determined  by  the  rank  of  the  observability  matrix  [13]: 

k  CA„  CA^  CA^]’'.  (13) 

From  (13),  the  system  is  observable  for  non-zero  angular  velocity.  This  result  supports 
the  well-known  drawback  of  the  sensorless  PMSM  -  its  inability  to  estimate  rotor  angle 
at  zero  and  low  angular  velocities  [14,  IS]. 

Given  the  ill-conditioning  of  the  rotor  position  observer  at  standstill,  special  provisions 
must  be  made  to  start  the  sensorless  PMSM  and  to  operate  the  system  at  extremely  low 
angular  velocity.  For  the  sensorless  drive  demonstrator  system,  the  observer  feedback  is 
disabled  at  speeds  below  7.5  r/s  (18  RPM),  which  represents  an  experimentally 
determined  low-speed  threshold. 
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2.3  Angular  Velocity  Estimation 


It  is  clear  from  (8)  that  rotor  velocity  is  a  required  parameter  for  the  implementation  of 
the  proposed  rotor  position  observer.  In  addition,  angular  velocity  feedback  is  required  to 
accurately  track  the  velocity  reference  in  a  speed  control  loop  used  in  most  propulsion 
systems.  Therefore,  in  the  absence  of  an  angular  velocity  sensor,  a  suitable  strategy  must 
be  developed  to  determine  this  parameter.  Two  strategies  are  described  in  the  following 
sections,  and  the  strengths  of  each  technique  are  combined  to  form  an  adaptive  velocity 
estimation  scheme. 


2.3.1  Quasi  Steady  State  Velocity  Estimation 

A  quasi  steady  state  estimate  for  magnitude  of  the  angular  velocity  [16]  is 


(14) 


This  result  suggests  the  use  of  the  input  and  output  vectors,  u  and  y,  to  estimate  the 
magnitude  of  the  angular  velocity.  The  direction  of  the  angular  velocity  estimate  at 
.«;ampling  interval  k  is  then  obtained  from 

dfg  =  sgn(^ik)  -  (f>{k  -  l))|<ye| ,  (15) 


where 


(^{k)  =  tan"* 


'  Ri^{k)-v„{k)' 

^Vp{k)-Rip{k)^ 


(16) 


The  strength  of  this  method  is  its  ability  to  determine  velocity  -  even  at  zero  and  very 
low  speeds.  The  weakness  is  its  dependency  on  PMSM  parameters. 


2.3.2  Time  Derivative  of  Estimated  Rotor  Angle 

An  alternative  to  the  method  given  in  (15)  is  to  use  the  time  derivative  of  the  estimated 
rotor  angle 


o>. 


(17) 


This  method  obviously  requires  a  stable  position  observer,  and  the  observer  pole 
locations  should  be  much  faster  than  the  angular  velocity.  Both  of  these  constraints  nwy 
be  enforced  by  proper  selection  of  the  system  eigenvalues  at  non-zero  angular  velocity. 
Subject  to  these  constraints,  the  estimated  angular  velocity  obtained  by  (17)  approaches 
the  actual  value  when  averaged  over  a  suflBcient  period  of  time,  as  shown  in  [18]. 

The  drawback  to  this  method  is  the  noise  that  results  from  the  differentiation  process  and 
its  inability  to  obtain  an  accurate  speed  estimate  at  startup,  where  the  rotor  position  is  not 
observable.  The  strength  of  this  method  is  the  accuracy  of  the  steady  state  speed 
estimates  throughout  the  medium  to  high  speed  range. 


5 


2.3.3  Adaptive  Velocity  Estimation 

The  two  proposed  velocity  estimation  methods  complement  each  other  nicely.  The  use  of 
(15)  yields  good  performance  at  very  low  speeds,  but  is  subject  to  parameter  uncertainty. 
The  use  of  (17),  however,  is  essentially  independent  of  the  parameter  estimation  errors  so 
long  as  the  system  poles  are  chosen  to  be  faster  than  the  angular  velocity  and  the  velocity 
is  sampled  over  a  sufficient  period  of  time.  Thus,  the  strengths  of  each  method  may  be 
combined  to  form  the  velocity  estimation  correction  scheme,  as  shown  in  Figure  1. 


Figure  1.  Velocity  estimation  correction  block. 


To  use  (17)  to  correct  for  velocity  estimation  errors  obtained  by  (15),  (14)  is  modified  to 
include  an  uncertainty  term  A,  which  represents  the  uncertainty  in  the  permanent  magnet 
flux  linkage: 


Initially,  the  error  term  A  is  set  to  zero.  However,  as  long  as  the  m^nitude  of  the 
estimated  velocity  is  above  the  low  speed  threshold  where  the  position  observer  is 
assumed  to  be  stable,  the  adaptive  velocity  estimation  process  is  activated. 

The  adaptation  process  compares  the  output  of  (18)  with  (17)  to  generate  the  velocity 
estimation  error.  A  PI  controller  operates  on  the  estimation  error,  adjusting  the  output  A 
to  correct  for  the  parameter  uncertainty.  The  time  constant  of  the  PI  controller  is  chosen 
to  be  slow,  since  the  temperature  effect  that  causes  the  parameter  drift  is  a  very  slow 
effect  as  well.  This  has  the  benefit  that  the  noise  generated  by  the  differentiation  of  the 
rotor  angle  estimate  is  smoothed  over  time.  Thus,  as  long  as  the  rotor  position  estimator 
remains  stable,  the  average  value  of  its  derivative  is  close  to  the  actual  rotor  speed. 

This  velocity  estimation  scheme  has  been  shown  through  experimentation  [17]  to 
produce  accurate  trackii^  of  the  velocity  reference,  despite  the  absence  of  the 
position/velocity  sensor.  This  is  of  critical  importance  in  undersea  vehicles  where  the 
rotor  angular  velocity  is  used  to  achieve  the  desired  forward  vehicle  speed.  An  overall 
block  diagram  of  the  proposed  position  observer,  including  the  angular  velocity 
estimation  block,  is  shown  in  Figure  2. 


6 


Figure  2.  Block  diagram  of  rotor  position  estimator. 


2.4  Rotor  Position  Observer  Pole  Placement 


From  the  system  model  given  in  (6),  the  open-loop  eigenvalues  of  the  plant 

(PMSM)  are 

^PMSM3i  ^PMSM4  -  ^J^e 

In  general,  the  location  of  the  desired  observer  eigenvalues  A  can  be  represented  by  their 
locations  on  the  real  (Re)  and  imaginary  (Im)  axis: 


A,  =  Re,  +  y  Imi  ;  =  A*  =  Re,  -  y  Im, 

Aj  =Re2  +  yimj  ;  =X\  =Re2- yini2 

The  system  eigenvalues  are  moved  from  their  open-loop  locations  (19)  to  the  desired 
locations  (20)  through  the  selection  of  the  observer  gain  matrix  G  in  (10).  It  is  the 
observer  eigenvalue  locations  that  determine  the  rate  of  convergence  of  the  angle 
estimate  error  towards  zero  and  the  stability  of  the  observer.  In  this  section,  the  relation 
between  the  observer  gain  matrix,  the  operating  condition  of  the  PMSM,  and  the  desired 
eigenvalues  is  presented. 

It  will  be  shown  later  that  it  is  advantageous  to  select  the  eigenvalues  of  the  observer  such 
that  their  magnitude  is  much  greater  than  the  PMSM  angular  velocity.  Therefore,  the 
observer  eigenvalues  are  placed  according  to  the  angular  speed  of  the  rotor  (Og : 
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(21) 


y^l  —  ^3  ^X-real  J^X-imag) 

^2  ~  ^4  ~  ^  1“  |®e  |(~  ^X-real  ~~  J^X-imag 

where  ^^-rea/  ^x-imag  define  the  response  and  damping  characteristics  of  the 

observer.  In  practice,  the  eigenvalues  of  the  position  observer  should  always  be  selected 
such  that  the  position  estimation  error  converges  at  a  fast  rate  with  respect  to  the  rotor 
angular  velocity.  Therefore,  in  the  experimental  system,  was  selected  so  that  the 

observer  poles  were  at  least  four  times  faster  than  the  rotor  velocity  -  that  is  ^  4 . 
The  imaginary  pole  placement  constant  was  selected  to  be  approximately  half  of 

^x-rea! »  which  was  found  to  produce  good  damping  characteristics.  Therefore,  from  (21), 
the  gain  scheduler  in  the  protot)'pe  sensorless  PMSM  drive  selects  the  observer  gain 
values  based  on  the  estimated  rotor  velocity. 

Each  time  a  new  velocity  estimate  is  available,  the  observer  gain  matrix  G  is  updated  by 
using  conventional  pole-placement  methods.  With  this  technique,  the  eigenvalues  of  the 
characteristic  equation 


|.II-(A„-GC)|  =  0 


(22) 


are  placed  in  the  left-half  plane  at  the  locations  specified  by  (21)  so  that  the  observer 
estimation  error  e  converges  to  zero  within  a  satisfactory  period  of  time.  In  [18],  a 
solution  to  (22)  is  developed,  and  the  observer  gain  matrix  G  required  to  place  the 
eigenvalues  at  the  location  specified  by  (21)  is  found  to  be: 


-R 

RX  1^2 

(OgT 

+  •^2) 


RX  1^2 
as 


-R 

X]X2  +fl>g  ) 


r 

X-,  +0)^ 


(OS 


Rix^  +  •<^2) 
r 


(23) 


Note  that  the  observer  gain  matrix  G  is  dependent  on  the  rotor  velocity,  the  desired 
eigenvalue  locations,  and  the  PMSM  stator  winding  parameters.  In  practice,  since 
velocity  is  estimated  in  a  sensorless  system,  the  gain  matrix  G  is  calculated  using 

instead  of  using  co^  as  indicated  in  (23). 
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2.5  Estimation  Error  Analysis 


The  error  dynamics  of  the  state  estimator  given  by  (10)  for  the  M-order  observer  shows 
that  asymptotic  state  reconstruction  is  achieved  if  the  eigenvalues  of  -GC  have 
negative  reals.  For  parameter  uncertainty,  the  system  matrix  of  the  plant  and  reference 
model  differ,  resulting  in  state  reconstruction  error.  In  this  section,  the  effects  of  such 
parameter  uncertainties  on  the  rotor  position  estimator  accuracy  are  investigated. 


To  quantify  the  effect  of  parameter  uncertainties  on  rotor  angle  estimation  accuracy,  (9) 
is  used  with  (1 1)  to  reveal  the  relation  between  the  rotor  angle  estimation  error  and  the 


error  associated  with  estimated  states  X3  and  : 


^  =  tan 


-1 


J^4-g4 

^3-^3 


=  tan 


-1 


^  ^pmfi  ^4 
y^pma 


(24) 


For  the  convenience  of  graphical  representation  (24)  can  be  expressed  as  a  space  phasor: 

~  ^4^  (■^pma  ~^3)  >  (25) 

where  A:  is  the  amplitude  of  the  rotating  phasor.  This  is  shown  graphically  in  Figure  3, 
which  clearly  shows  how  state  estimation  errors  ^3  and  result  in  a  difference  between 

A 

the  actual  and  estimated  rotor  angles  -  shown  as  ^  and  ^ ,  respectively. 

In  the  event  of  zero  estimation  error,  (24)  shows  that  the  estimation  error  is  equal  to  the 
actual  angle  since 


^  =  tan  ‘ 


^  yl3/2A„sm0^ 

COS0 J 


^9  . 


(26) 


real 


Figure  3.  Effect  of  state  estimation  error  on  angle  estimation. 
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In  this  case,  the  space  phasor  defined  by  (25)  is  defined  in  steady  state  by  a  circle  of 
radius  ^  centered  on  the  origin.  In  the  case  of  nonzero  state  estimation  error, 

however,  the  trajectory  traced  by  the  space  phasor  of  (25)  is  no  longer  defined  by  a  circle 
centered  at  the  origin.  The  trajectory  of  this  space  phasor  is  a  convenient  graphical  tool 
with  which  to  characterize  the  rotor  angle  estimation  error,  and  will  be  used  in  the 
following  sections  with  the  goal  of  characterizing  the  effects  of  parameter  uncertainties 
on  the  rotor  angle  estimation  accuracy. 

2.5.1  Effect  of  Angular  Velocity  Estimation  on  Rotor  Angle  Estimation  Accuracy 

For  a  shaft  sensorless  PMSM  drive,  rotor  velocity  must  be  estimated  as  well  as  the  rotor 
position.  It  follows  that  the  slowly  varying  velocity  and  its  estimated  value  Sg  are 

related  by 


<i)g  =Sg+ACDg  , 


(27) 


where  is  the  velocity  estimation  error.  From  the  state  model  of  (6),  where  A^  is 
velocity  dependent,  the  plant  behavior  can  be  expressed  with  respect  to  the  estimated 
velocity  and  the  corresponding  velocity  estimation  error: 


x  =  A^x  +  A^^x  +  Bu  , 


(28) 


where 


-  r 

0 

T 

0 

'0 

0 

0 

0 

0 

-  r 

0 

r 

0 

0 

0 

0 

> 

li 

0 

0 

0 

5 

0 

0 

0 

-  Ad) g 

0 

0 

0 

0 

0 

ACDg 

0 

(29) 


Similarly,  the  reference  model  will  behave  accordit^  to 


i  =  A^x  +  Bu  +  GC(x-x)  . 


(30) 


From  the  state  estimation  error  given  in  (11), 

e  =  x-i.  (31) 

so  that  from  (28)  and  (30),  the  error  dynamics  with  uncertainty  in  the  angular  velocity 
become 


e  =  (A^  -GC)e  +  AA^x  . 


(32) 
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In  order  for  the  errors  in  (32)  go  to  zero  asymptotically,  independent  of  x,  the  coefBcient 
of  X  would  have  to  be  zero  and  -GC  must  be  the  dynamics  matrix  of  a  stable 

system.  When  these  conditions  are  met,  the  error  dynamics  are  reduced  to  the 
homogeneous  equation,  and  the  error  approaches  zero  in  steady  state.  In  practice,  the 
requirement  for  zero  angular  velocity  estimation  error  is  a  difficult  objective  and 
therefore  an  analysis  of  velocity  uncertainty  on  angle  estimation  accuracy  is  required. 

In  order  to  quantify  the  effects  of  the  velocity  estimation  error  on  rotor  position 
estimation  accuracy,  (32)  is  solved  to  reveal  the  steady  state  estimation  error 
corresponding  to  estimated  states  Xj  and  : 


f 

2  ^ 

sin(^(0  +  <^) 

Sg  i 

-  cos(^(r)  +  p) 

(33) 


where 


r  col-^Ax  ^ 

(34) 


The  computational  details  of  this  result  are  included  in  [18].  Now,  the  estimated  states 
with  velocity  estimation  uncertainty  can  be  written  in  terms  of  the  estimated  error: 


COS^(0  = 


(Oe  ) 

f  A  ^ 

sin  0(t)  +  —^Kx  cos(^(/)  +  ip) 
(o.  ; 


(35) 


which  through  some  trigonometric  manipulations  can  be  expressed  by; 


'  A 

^3 

X4 


[3,  .icos(0(t)+<^y 

h  ”  lsm(0(t)+Sl 


M  = 


tr  X 


\ 


+ 1- 


bkOi. 


CO, 


5  =  tan 


-1 


CO. 


\ 

Kx  sin^ 

e  ) 


Kx  cosp 


l-^Kxsmp 


(36) 
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Substitution  of  (36)  into  (9)  reveals  the  effect  of  the  velocity  estimation  error  on  position 
estimation  accuracy; 


S{t)  =  tan‘‘ 


^  .J3/2A„M  sin(^(f)  +  S)  ^ 


=  0(t)  +  S 


(37) 


where  S  is  the  position  estimation  error  that  results  from  angular  velocity  uncertainty. 
This  result  implies  that  for  a  given  observer  gain  G,  the  position  estimation  accuracy  is 
bound  by  the  accuracy  of  the  velocity  estimation  Furthermore,  it  demonstrates  that 
selection  of  the  system  eigenvalues  such  that  |®e|  (i'®->  choose  and 

ki-imag  ®  ^  much  greater  than  one)  reduces  the  effects  of  velocity  estunation 

error  by  forcing  ^  -  tc/I  ,  so  that  — >  0 . 

Figure  4  shows  the  locus  of  the  estimated  state  and  angle  estimation  error  resulting  from 
a  ±  10%  velocity  estimation  error  and  positive  angular  velocity.  To  accentuate  the  angle 
estimation  error,  the  closed-loop  eigenvalue  reals  were  approximately  equal  to  the 
angular  velocity.  In  Figure  4,  is  the  angle  between  a  line  drawn  from  the  origin  to  the 
point  on  the  locus  marked  ^  =  0  and  the  Xj  axis  on  the  locus  plots,  or  the  vertical 

distance  0-0  in  the  plots  showing  actual  and  estimated  angles.  The  estunated  angle 
leads  and  lags  the  actual  ai^le  for  negative  and  positive  velocity  estimation  errors, 
respectively* 

The  effect  of  moving  the  pole  locations  further  into  the  left  half  plane  is  shown  in  Figure 
5.  Increasing  the  magnitude  of  the  conjugate  pair  Aj,  ^  drives  the  angle  <ff  toward 

and  5  toward  zero.  Figure  5  represents  the  same  operating  condition  as  the  second 
column  in  Figure  4,  but  in  conyarison  has  much  lower  angle  estimation  error  due  to  the 
pole  locations. 
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Negative  Vel.  Est  Error  Positive  Vel.  Est  Error 


Figure  4.  Angle  estimation  error  due  to  ±10%  velocity  error,  slow  observer  poles. 
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Figure  5.  Angle  estimation  with  10  percent  velocity  error,  fest  observer  poles. 


2.5.2  Effect  of  Stator  Resistance  Uncertainty  on  Angie  Estimation  Accuracy 

In  the  PMSM,  the  most  temperature  sensitive  parameter  is  the  stator  winding  resistance, 
which  varies  linearly  over  the  operating  temperature  of  the  machine  [19].  It  is  therefore 
reasonable  to  assume  that  the  angle  estimation  process  will  be  subject  to  uncertainty  in 
this  parameter,  and  a  study  of  its  effect  is  crucial. 
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The  relation  between  the  actual  stator  resistance  and  the  nominal  value  used  by  the 
reference  model  can  be  represented  by 

R  =  R-^-AR  ,  (38) 


where  R  is  the  stator  resistance  value  used  by  the  PMSM  internal  model,  R  is  the  actual 
PMSM  stator  resistance,  and  M  represents  the  modeling  error.  By  using  (38)  and 
assuming  that  all  other  parameters  are  accurately  known,  the  PMSM  system  matrix  given 
in  (8)  can  be  expanded  to  reflect  the  resistance  uncertainty: 


^<0  ^  ^AR 


R 

R 

AR 

AR 

0 

0 

0 

— 

0 

'X 

Ls 
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R 

0 

A 
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0 

AR 
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Ls 

0 
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0 
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0 

0 

0 

0 

0 

(39) 


Thus,  the  PMSM  and  reference  model  behave  accordii^  to 


x  =  A^x  +  A^x  +  Bu 
X  =  A^i  +  Bu  +  GC(x  -  i) 


(40) 


respectively,  so  the  error  dynamics  defined  in  (31)  become 

e  =  x-x  =  (A<„ -Gc)e+A^x  . 


(41) 


The  solution  to  (41)  fi-om  [18]  reveals  the  steady  state  estimation  errors  introduced  in  the 
states  Xj  and  JC4  due  to  resistance  uncertainty: 


AR 

-Ap  sm(0  +  s)+J 

3,  •  / 

—A-  sm  ofj 
2  ” 

0>eLs 

ApCOs(ff  +  s)-J^ 

~A„  cosOgt 

(42) 


In  (42),  the  angle  s  represents  the  angle  between  the  resultant  stator  flux  and  the  flux 
generated  by  the  rotor  permanent  magnets  alone.  For  the  lightly  loaded  or  low- 
inductance  PMSM  such  as  a  slotless  design,  this  angle  will  be  approximately  zero.  For  a 
smaller  airgap  machine,  the  angle  e  will  increase  with  stator  current. 
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From  the  state  estimation  error  in  (42),  the  effects  on  the  rotor  angle  estimation  accuracy 
due  to  resistance  uncertainty  are  determined  from  (9),  (11),  and  (42): 


9{t)  =  tan 


-1 


■^4  ~  ^4,M 
V  ^3  ~  J 


=  tan" 


^  [3  ^ 

J-/l„sin^(0-e4A/j 

^^cos^(r)-e3^ 


(43) 


Substitution  of  (42)  into  (43),  and  the  trigonometric  manipulations  detailed  in  Appendix 
B  of  [18]  yields  the  angle  estimation  error  induced  by  resistance  uncertainty: 


9^  =tan-' 


^  sin(^  + 
cos(9  +  <^/} ) 


=  9  +  S, 


=tan  ‘ 


AR 


(4sin^) 


cos^l 


(o. 


(44) 


where  4  represents  the  angle  of  the  stator  current  with  respect  to  the  generated  emf. 
From  the  result  of  (44),  the  effect  of  stator  resistance  uncertainty  is  to  introduce  a  steady 
rotor  angle  estimation  error  of  Sn ,  whose  magnitude  is  proportional  to  the  resistance 
uncertainty  and  inversely  proportional  to  angular  velocity.  The  error  is  shown  in  Figure 
6,  where  is  the  angle  between  the  x-^  axis  and  a  line  drawn  from  the  origin  to  the  point 
on  the  locus  marked  ^  =  0 .  In  the  simulation  shown  in  Figure  6,  the  nominal  resistance 
is  1.35  while  the  actual  stator  resistance  due  to  heating  is  2.35  O .  Angular  velocity  is 
100  r/s,  and  the  machine  is  assumed  to  be  operating  such  that  an  angle  of  ^  =  -;r/10  is 
present  between  the  back  emf  and  the  stator  currents. 

In  the  above  error  analysis,  it  was  assumed  that  the  resistance  uncertainty  is  balanced.  In 
the  event  of  unbalanced  resistance  error,  the  locus  shown  in  Figure  6  is  deflected  into  an 
oval  shape.  The  non-circular  locus  of  the  estimated  states  results  in  an  additional  angle 
estimation  error  at  twice  the  fundamental  frequency  in  addition  to  the  constant  ai^le  error 
term  5 ^ . 
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Figure  6.  Resistance  uncertainty  effect  on  angle  estimation:  A/?  -  IQ ,  heavily  loaded. 

In  the  expression  for  the  numerator  is  seen  to  be  dependent  on  the  magnitude  of  the 
stator  current.  In  the  extreme  case  of  zero  stator  current,  the  numerator  goes  to  zero  and 
no  angle  estimation  errors  resuk  from  resistance  uncertainty.  This  is  an  easily  justifiable 
resuk,  since  angle  estimation  of  an  open  circuk  stator  PMSM  does  not  require  knowledge 
of  the  stator  resistance.  It  should  also  be  noted  that  if  the  stator  current  is  in  phase  with 
the  back  emf  (i.e.  ^  =  0),  uncertainty  in  the  stator  resistance  generates  no  error  in  the 

estimated  rotor  angle. 

2.5.3  Effect  of  Stator  Inductance  Uncertainty  on  Angle  Estimation  Accuracy 

In  the  event  that  the  stator  inductance  is  not  known  accurately,  an  angle  e^imation  error 
will  resuk.  The  characteristics  of  the  angle  estimation  error  due  to  inductive  uncertainty 
are  now  analyzed  by  representing  the  actual  stator  inductance  as 

L,=Z,  +  AI,  (45) 

where  L^,  I,,  and  AL  are  the  actual  value,  estimated  value,  and  uncertainty  of  the  stator 
inductance,  i,  =1,  Given  the  inductance  uncertainty  defined  in  (45),  the 

PMSM  dynamics  may  be  described  by 

x  =  A^x  +  A^x  +  Bu  ,  (46) 

where  the  actual  system  matrix  given  in  (8)  is  expanded  to  include  the  inductance 
uncertainty  term: 
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The  system  matrix  used  by  the  PMSM  reference  model  and  the  corresponding  uncertainty 
in  the  system  matrix  due  to  inductance  error  matrix  in  this  case  are  described  by 


respectively. 

In  addition,  the  inductance  error  also  appears  in  the  output  matrix  C  defined  in  (8).  To 
include  this  effect,  the  PMSM  output  matrix  defined  in  (8)  can  be  rewritten  to  account  for 
the  uncertainty  in  the  stator  inductance: 


C  =  C  +  C^  - 
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(49) 


where  the  output  matrix  used  by  the  reference  model  is 
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and  the  error  in  the  output  matrix  created  due  to  the  uncertainty  in  the  inductance  is 


Ca,= 


AL 

4(4 +Ai) 
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(51) 


The  PMSM  reference  model  dynamics,  in  the  case  of  inductance  uncertainty,  are 

i  =  A^i  +  Bu  +  G{(c  +  C^)x-Ci},  (52) 

and  finally,  the  error  dynamics  are  determined  from  (46)  and  (52)  are 

e  =  x-i  =  (A,-Gc)c+(A^-GC^)x  .  (53) 

From  Appendix  B  of  [18],  the  steady  state  solution  to  (53)  is 
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The  effects  of  the  state  estimation  error  caused  by  inductance  modeling  uncertainty  may 
be  related  to  angle  estimation  accuracy  by  (9)  and  (11).  In  Appendix  B  of  [18],  it  is 
shown  that  the  estimated  state  variables  can  be  written  as 

Xj  =  J— cos^(0 - ^3  =  cos{0  +  Si) 

_  ,  (55) 

^4  =  sin  «(/)-«, 


18 


0.2 


0 


-0.1 

-0.15 


-0.2*“ 

0 


0.02  0.04  0.06  0.08 

tjme(s) 


Figure  7.  Stator  inductance  uncertainty  effect  on  rotor  angle  estimation. 


where 
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=  ALI„cos4  ;d  =  -  ALI„  sin^ 
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(56) 


From  this  result,  inductance  error  introduces  a  constant  angle  estimation  error  Si ,  whose 
amplitude  is  proportional  to  the  magnitude  of  the  inductance  uncertainty,  the  stator 
current  an^litude,  and  the  angle  ^  of  the  stator  current  with  respect  to  the  generated  emf. 
The  result  of  inductance  uncertainty  is  shown  in  the  locus  of  Figure  7,  where  the  offset 
Si  is  illustrated  by  the  ai^le  between  the  x,  axis  and  a  line  drawn  from  the  origin  to  the 
point  marked  as  ^  =  0 .  In  the  event  of  unbalanced  inductance  xmcertainty,  the  locus  in 
Figure  7  is  distorted  such  that  an  additional  double  frequency  term  is  added  to  the 
constant  estimation  error  Si . 


2.5.4  Effect  of  PM  FIox  Linkage  Uncertainty  on  Angle  Estimation  Accuracy 

Although  the  PM  flux  linkage  constant  is  not  a  parameter  in  the  PMSM  reference  model 
used  for  rotor  angle  estimation,  it  is  required  to  estimate  angular  velocity  when  the  slow 
mode  approach  is  used.  Therefore,  when  this  method  of  velocity  estimation  is  active,  the 
effect  of  the  permanent  magnet  flux  linkage  xmcertainty  can  be  considered  to  be  the  same 
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as  a  rotor  velocity  estimation  error,  since  uncertainty  in  this  parameter  translates  to 
uncertainty  in  the  angular  velocity. 


2.5.5  Effect  of  Stator  Voltage  and  Current  Measurement  Errors 

The  implementation  of  the  rotor  position  observer  requires  the  acquisition  of  the  stator 
currents  and  vohages.  The  use  of  these  signals  in  the  rotor  position  estimator  will  result 
in  estimation  inaccuracy  when  the  sensed  quantities  contain  scaling  or  offset  error. 
Furthermore,  the  inaccuracies  may  be  unequal  or  unbalanced  from  one  phase  to  another. 
The  effects  of  these  measurement  inaccuracies  are  examined  in  this  section. 


2.5.5.1  Voltage  Measurement  Error 

In  the  case  of  voltage  sensing  error,  the  actual  and  acquired  signals  are  related  by 

^c^.meas  ~  .actual  afi  *  (^”7) 


In  this  case,  the  PMSM  dynamics  can  be  expanded  to  include  the  measurement  error: 

X  =  A^x  +  b(v^  -  AV^ )  ,  (58) 


and  the  reference  model  becomes 

i  =  +  +GC(x-x)  .  (59) 


Therefore,  the  error  dynamics  are 

e  =  (A,-GC)e-BAV^  .  (60) 


From  [18]  the  steady  state  solrition  of  (60)  for  balanced  voltage  scaling  errors  and 
unbala^ed  voltage  offset  errors  is: 
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respectively,  where  AV,  Av^,  and  Av^  are  the  magnitudes  of  the  balanced  voltage 
scaling  error,  and  the  offset  error  in  the  fictitious  two  phase  windings. 

For  a  balanced  scaling  error  on  the  voltage  sensor,  the  solution  for  the  resulting  angle 
estimation  error  is  determined  in  [18],  which  gives: 


A 

e 


=  tan 


=  0-^5, 


v,scaie 


(63) 


This  result  clearly  shows  that  the  resulting  angle  estimation  error  is  a  constant  ai^le  for  a 
balanced  sensor  scaling  inaccuracy.  Under  these  conditions,  the  circular  locus  of  the 
estimated  states  is  rotated  by  the  ai^le  1*^®  the  locus  shown  previously  for 

stator  resistance  and  inductance  uncertainties. 


In  the  case  that  the  scaling  error  is  not  balanced  between  the  phases,  the  estimated  state 
locus  becomes  non-circular,  and  an  angle  estimation  error  of  double  the  fundamental 
rotational  frequency  is  superimposed  upon  the  constant  estimation  error  angle  . 
The  estimated  state  locus,  estimated  states,  and  angle  estimation  error  are  each  shown  in 
Figure  8  for  an  unbalanced  voltage  scaling  error. 

In  the  case  of  voltage  offset  error,  examination  of  (62)  reveals  that  a  constant  error  that  is 
not  a  function  of  the  rotor  angle  is  introduced  into  the  estimated  state  variables.  This 
resuhs  in  an  estimated  state  locus  that  is  not  centered  at  the  origin.  Furthermore,  if  the 
voltage  measurement  offsets  are  not  balanced,  the  locus  becomes  non-circular,  as  shown 
in  Figure  9. 

The  translated  locus  in  Figure  9  generates  an  angle  estimation  error  at  the  fundamental 
frequency,  while  a  non-circular  locus  creates  an  additional  error  term  at  twice  the 
fundamental  frequency.  The  peak  magnitude  of  the  angle  estimation  error  in  the  case  of 
unbalanced  voltage  offset  is 
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pyVoffset 
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^3ss,Voff  ^Ass,Voff 


(64) 
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Thus,  from  (62)  and  (64)  it  can  be  seen  that  selecting  the  magnitude  of  the  closed-loop 
eigenvalue  reals  to  be  much  greater  than  the  angular  velocity  helps  to  reduce  the  effect  of 
voltage  measurement  offeets  on  the  rotor  position  estimation  accuracy. 


time  (s) 


Figure  8.  Effect  of  unbalanced  voltage  sensor  scaling  error. 


Figure  9.  Effect  of  voltage  measurement  offset  of  angle  estimation. 
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2.5.5.2  Current  Measurement  Errors 

In  the  case  of  current  measurement  inaccuracies,  the  actual  and  measured  signals  are 
related  by 


^afi,aieas  ~  ^afi,actuai  ^aft  » 

Where  >  Ja^.actuai^  ^nd  represent  the  measured  and  actual  stator  currents,  and 

the  error  associated  with  the  measurement,  respectively.  With  the  current  measurement 
relations  defined,  the  output  equation  (7)  is  modified  to  account  for  the  current 
measurement  inaccuracy: 


y  =  cx  =  v, 


actual 


—  A/^ 

^a,meas  a 
^ fi.msta  ~ 


(66) 


The  reference  model  may  be  expanded  to  include  the  current  measurement  uncertainty: 

i  =  Aa,i  +  Bu  +  G((l^,j^ +AlQjg)-i(^)  .  (67) 

Thus,  the  error  dyiramlcs  are 


e  =  (A^-GC)e-GAI^  . 


(68) 


From  Appendix  B  of  [18],  the  steady  state  solutions  to  (68)  for  current  scaling  and  offset 
errors  are: 
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(70) 


The  terms  A/ ,  A/^ ,  and  Mp  represent  the  magnitude  of  the  balanced  current  scaling 
error,  and  current  offsets  in  the  a  and  axes,  respectively. 
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The  effect  of  current  sensor  scaling  errors  on  the  angle  estimation  accuracy  is  computed 
in  Appendix  B  of  [18].  For  the  balanced  current  scaling  error,  the  introduced  ai^le 
estimation  error  is 


^  =  tan  ^ 
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(71) 


Therefore,  a  balanced  current  scaling  error  introduces  a  constant  angle  estimation  error  of 
Si  ,  whose  magnitude  is  dependent  on  the  magnitude  of  the  scaling  error.  From  (71), 

the  angle  estimation  error  is  more  prone  to  current  scaling  errors  at  low  speeds.  Similar 
to  the  voltage  measurement  scaling  error,  any  unbalance  in  the  current  scaling  error 
vector  will  generate  a  double  frequency  ar^le  estimation  error  term  that  is  superimposed 
onto  the  angle  Sf  ^aie. 

From  (71),  the  effect  of  current  offset  error  is  a  constant  term  added  to  the  estimated 
states.  Thus,  current  offset  is  identical  in  nature  to  the  voltage  offset  error  described  in 
the  previous  section  with  the  exception  that  the  amplitude  of  the  angle  estimation  error  in 
the  cxirrent  offeet  case  is 
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Assjoff 


(72) 


2.6  Eigenvalue  Selection  Strategy 

The  observer  poles  should  be  selected  according  to  the  system  requirements,  which  are 
typically  r^id  response  and  stable  estimation.  As  the  real  parts  of  the  eigenvalues 
become  larger  negative  values,  the  estimation  error  converges  rapidly.  Choosing  the 
eigenvalues  too  aggressively,  however,  results  in  the  an^)lification  of  system  noise  and 
reduced  or  even  unstable  performance. 

A  general  a^orhhm  for  observer  pole  placement  was  presented  in  (21),  where  it  is 
suggested  that  the  observer  poles  be  placed  based  on  the  rotor  velocity.  Furthermore,  the 
results  of  the  angle  estimation  error  study  indicate  that  choosing  observer  poles  much 
faster  than  the  angular  velocity  of  the  PMSM  reduces  the  effects  of  angular  velocity 
estimation  uncertainty  as  well  as  voltage  and  current  measurement  errors.  Velocity 
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estimation  xising  the  time  derivative  of  the  estimated  rotor  angle  is  also  improved  by 
choosing  feist  observer  eigenvalues  with  respect  to  angular  velocity,  as  shown  in  [18]. 

In  the  sensorless  drive  technology  demonstrator  system,  the  eigenvalues  of  the  position 
observer  were  scheduled  according  to  the  angular  velocity  as  indicated  in  (21).  Through 
experimentation,  it  was  found  that  observer  poles  at  least  4  times  fester  than  the  angular 
velocity  produced  excellent  results.  Thus  in  (21),  the  pole  placement  was  performed  with 
^x-reai  ^  4  and  The  real  part  is  chosen  for  stability  purposes  while 

the  imaginary  part  was  foimd  through  experimentation  to  produce  desirable  damping 
characteristics. 


3  SENSORLESS  DRIVE  IMPLEMENTATION 

In  this  chapter,  the  rotor  position  estimation  algorithm  outlined  in  Chapter  2  is 
implemented  on  a  digital  signal  processor  (DSP).  The  processing  hardware,  associated 
power  electronics,  and  sensors  used  to  realize  the  sensorless  drive  technology 
demonstrator  are  presented  and  implementation  details  are  discussed. 


3.1  Control  Implementation 

Electric  drives  are  used  to  control  the  torque  and  speed  of  a  motor  in  a  closed  loop 
system.  A  resolver  or  encoder  attached  to  the  shaft  of  the  motor  is  normally  used  to 
provide  the  feedback  required  to  carry  out  the  torque  and  velocity  control  loops.  In  the 
sensorless  drive  technology  demonstrator,  however,  the  rotor  position  and  velocity  are 
supplied  by  the  rotor  angle  and  velocity  estimator,  respectively. 

The  arrangement  of  the  sensorless  controUer  is  shown  in  Figure  10,  where  cascaded 
torque  and  speed  loops  are  shown.  The  actual  drive  is  shown  in  Figure  11.  The  outer 
velocity  loop  generates  the  torque  reference  required  to  drive  the  speed  of  the  PMSM 
shaft  toward  its  commanded  value.  From  the  torque  command  and  knowledge  of  the 
rotor  angle,  balanced,  three  phase  current  commands  that  are  required  to  produce  the 
commanded  torque  are  generated.  The  PWM  current  controller  then  applies  the  required 
voltage,  via  the  inverter,  to  the  PMSM  terminals  in  order  to  enforce  the  current 
references. 
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3.1.1  Current  Controller 

The  innennost  torque  control  loop  is  clearly  the  most  important  loop  in  the  cascaded 
control  structure  shown  in  Figure  10.  In  order  to  obtain  optimum  performance  of  the 
entire  system,  including  the  outer  velocity  loop,  requires  that  the  developed  instantaneous 
torque  is  equal  to  the  commanded  torque.  In  the  sensorless  drive  system,  the  torque 
controller  is  implemented  through  a  current  regulator,  since  the  developed  torque  is 
directly  proportional  to  the  quadrature  axis  current  in  the  PMSM  [20].  The  current 
regulator  consists  of  a  voltage  source  PWM  inverter,  where  the  sinusoidal  PWM  method 
[21]  is  used. 

In  the  current  regulator  of  the  sensorless  electric  drive,  a  PI  controller  in  the  rotating 
reference  frame  has  been  adopted  [22,  23].  Current  control  is  greatly  simplified  in  the 
rotating  reference  frame,  making  a  three-phase  machine  equivalent  to  analyzing  a  DC 
motor.  The  rotating  reference  frame  current  controller  operated  on  DC  quantities,  and 
can  in  theory  regulate  the  currents  with  zero  steady-state  error. 

The  current  controller  implemented  in  the  sensorless  drive  is  based  on  the  technique 
presented  in  [22].  The  machine  model  in  the  synchronous  reference  frame  is  given  as 
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(73) 


In  (73),  vj  g  and  ^  represent  the  motor  inputs  and  outputs,  respectively  in  the  direct 

and  quadrature  axes.  In  the  synchronous  reference  frame,  the  induced  emf  in 

the  quadrature  axis  current  is  a  slowly  varying  disturbance  proportional  to  angul^ 
velocity.  In  addition,  (73)  shows  a  cross  coupling  of  the  direct  and  quadrature  axis 
currents  into  the  quadrature  and  direct  axis  voltage  expressions,  respectively.  If  the 
disturbance  and  cross  coupling  terms  are  decoupled  from  (73),  the  dependence  on  angular 
velocity  disappears  and  the  stator  voltage  equation  becomes  that  of  an  RL  (resistance  and 
inductance)  load,  enabling  simple,  fast,  and  accurate  current  regulation. 

Through  feedback  linearization,  a  decoupled  and  linear  system  can  be  obtained  by 
choosing 


(74) 


where  and  are  the  modified  output  voltages  of  the  current  regulator,  and  v,  are 
the  actual  voltages  applied  to  the  PMSM,  and  Vj  and  v“  are  the  decoupling  terms  : 
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In  (75),  the  decoupling  vector  is  used  to  overcome  the  rotor  velocity-dependent  cross 
coupling  and  motor  back  emf  terms.  Note  that  in  a  sensorless  drive,  the  calculation  of  the 
cross-coupling  term  must  use  the  estimated  rotor  velocity  ,  With  the  cross  coupling 
vector  defined,  the  system  may  be  described  by 
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to  which  linear  control  design  techniques  can  be  applied.  In  the  prototype  sensorless 
drive,  a  PI  controller  is  utilized  such  that 
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The  proportional  and  integral  gains,  and  K, ,  are  selected  to  be  equal  in  the  d  and  q 

axis  current  controllers,  and  are  determined  experimentally  to  give  a  closed  loop  current 
controller  of  sufficient  bandwidth  -  approximately  1  kHz  in  the  demonstrator  system. 


A  functional  block  diagram  of  the  current  control  used  on  the  prototype  sensorless  drive 
is  shown  in  Figure  12.  In  a  velocity  control  application,  an  outer  velocity  controller 
generates  the  current  or  torque  command  Tq* ,  which  is  translated  to  d-q  current 
commands  by  the  reference  current  generator.  The  direct  axis  current  reference  is 
normally  maintained  at  zero,  since  does  not  produce  torque  in  a  non-salient  PMSM.  That 
is,  the  quadrature  axis  current  command  /*  represents  the  primary  torque-producing 
component. 

In  Figure  12,  the  feedback  currents  are  transferred  into  the  rotating  reference  fi'ame 
through  Park  and  Clarke  transforms  and  the  rotor  angle,  which  is  generated  by  the  rotor 
position  estimator.  Con^arison  of  the  current  reference  to  the  feedback  yields  the  current 
error,  which  is  operated  on  by  the  PI  controller  to  generate  the  voltage  references  in  the 
rotating  reference  flame.  Finally,  the  voltage  reference  is  transformed  to  voltage 
commands  in  the  stationary  flame,  which  are  operated  on  by  the  PWM  controller  to 
generate  the  inverter  switching  control  signals.  Only  two  current  values  are  measured  by 
the  DSP  since  the  third  phase  current  is  not  an  independent  variable  for  a  Y-connected 
PMSM. 
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In  order  to  demonstrate  the  effectiveness  of  the  current  control  scheme,  some 
experimental  results  are  now  presented.  In  the  experiments,  the  proposed  current 
controller  is  implemented  on  a  digital  signal  processor  in  a  closed  loop  system  with  a  DC 
link  voltage  of 650  volts. 

In  Figure  13,  the  current  control  performance  in  steady  state  is  shown  for  an  operating 
speed  of  300  RPM,  and  load  torque  of  approximately  200  foot-pounds.  Figure  14 
demonstrates  the  transient  performance  during  a  step  change  in  the  quadrature-axis 
current  command,  which  occurs  at  approximately  200  i^liseconds  when  a  speed  reversal 
is  requested.  This  results  show  that  the  current  controller  exhibits  very  fast  control 
response  without  overshoot. 

For  the  startup  transient  shown  in  Figure  15,  the  quadrature  axis  current  is  required  to 
increase  with  the  speed  due  to  the  viscous  loading  on  the  system.  In  the  figure,  it  can  be 
seen  that  any  errors  in  the  estimated  rotor  angle  can  produce  unexpected  results,  since 
accurate  rotor  angle  feedback  is  required  to  correctly  control  the  torque  via  the  quadrature 
axis  current 
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CunwTt  Control  During  Startup  -  Viscous  Load 


Time  (s) 

Figure  15.  Current  control  during  startup  transient. 


3.1.2  Velocity  Controller 

In  order  to  implement  velocity  control,  a  proportional-integral  (PI)  speed  control  ^  been 
implemented.  Since  this  is  a  sensorless  drive,  the  velocity  controller  utilizes  the 
estimated  angular  velocity  (described  in  Chapter  2)  as  the  feedback.  The  PI  controUer 
operates  on  the  speed  error  to  generate  the  torque  command,  as  depicted  in  Figure  10. 

Steady  state  speed  control  performance  in  the  sensorless  mode  is  demonstrated  in  the 
experiment  Ulustrated  in  Figure  16,  where  the  actual  velocity  is  obtained  from  an  encoder 
signal  tettqwrarily  attached  to  the  shaft  of  the  motor  under  test.  In  this  experiment,  the 
commanded  velocity  is  -300  RPM,  the  mean  value  of  the  estimated  velocity  is  -300.1016 
RPM,  and  the  mean  value  of  the  actual  rpm  is  —300.0487  RPM.  The  resuhs  shown  in 
Figure  16  reveal  that  any  error  associated  with  the  estimated  angular  velocity  generates  a 
corresponding  ripple  in  the  actual  speed.  A  2nd  order  filter  with  a  cutoff  frequency  of  50 
r/s  conditions  the  estimated  angular  velocity.  While  a  lower  cutoff  frequency  would 
reduce  the  an^litude  of  the  ripple  on  the  estimated  velocity,  it  would  also  have  an 
adverse  effect  on  transient  performance. 

Figure  17  demonstrates  that  during  a  speed  reversal  transient,  filtering  of  the  estimated 
angular  velocity  results  in  a  delay  between  the  actual  and  estimated  angular  velocity. 
This  delay  is  seen  to  generate  difficulties  in  angle  estimation  -  especially  when  the 
estimated  angular  velocity  is  opposite  in  polarity  to  the  actual  estimated  an^lar  velocity. 
In  general,  if  the  estimated  and  actual  angular  velocity  are  of  opposite  sign  only  when 
observer  feedback  is  disabled  (i.e.  when  operating  below  the  low  speed  threshold),  the 
estimated  angle  error  will  recover  from  the  error  in  angular  velocity.  This  is  shown  in  the 
figure  where  the  angle  estimation  error  converges  toward  zero  when  the  angular  velocity 
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exceeds  the  low  speed  threshold  of  19  RPM  (at  approximately  1  second).  However,  if 
the  estimated  and  actual  angular  velocities  are  of  opposite  polarity  while  operating  above 
the  low-speed  threshold  level,  the  observer  becomes  unstable.  This  is  clearly  an 
undesirable  situation,  and  illustrates  the  considerations  in  selecting  the  cutoff  frequency 
of  the  estimated  angular  velocity  filter.  The  selection  of  this  cutoff  frequency  represents 
a  tradeoff  between  steady-state  and  transient  velocity  control  performance. 


T1me(s) 

Figure  16.  Sensorless  velocity  control  -  steady  state  300  rpm. 


Figure  17.  Sensorless  velocity  control  -  speed  reversal. 
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3.2  Hardware  and  Algorithm  Implementation 

In  this  section,  the  implementation  of  the  rotor  position  estimation  strategy  in  a  closed- 
loop  drive  system  is  presented.  Although  many  sensorless  control  algorithms  have  l^n 
proposed  in  literature,  very  few  have  actually  demonstrated  the  drive  performance  in  a 
closed-loop  system.  TTierefore,  some  issues  and  strategies  related  to  the  unplementation 
are  described. 


3.2.1  Computing  Hardware 

The  architecture  for  the  computing  hardware  of  the  prototype  sensorless  electric  drive 
consists  of  a  floating  point  digital  signal  processor  (DSP)  and  a  fixed  point  motor  control 
DSP.  A  block  diagram  of  the  arrangement  is  shown  in  Figure  1 8,  and  photo^phs  of  the 
controller  board  are  included  in  Figure  19.  In  the  demonstrator  sensorless  drive  system,  a 
personal  computer  acts  as  the  user  interfece,  providing  a  graphical  user  interfece  for 
operating  the  drive,  displaying  the  drive  performance  status,  and  storing  relevant 
operational  data.  The  developed  confuting  hardware  has  been  designed  to  perform 
sensorless  control  for  a  specific  Integral^  Motor/Propulsor  drive  application  [17].  Both 
DSPs  are  bootable  from  a  programmable,  erasable  memory  device  (EPROM). 

The  floating  point  DSP  (Analog  Devices  model  EZ-Kit  Lite)  performs  the  rotor  position 
and  velocity  estimation  at  100  ps  and  2  ms  intervals,  respectively.  The  fixed  point  motor 
control  DSP  (Analog  Devices  ADMC401)  is  used  to  implement  the  current  control 
scheme  previously  described,  sample  the  stator  voltages  and  currents,  and  provide  the 
PWM  signals  to  the  inverter.  The  PWM  frequency  for  the  sensorless  drive  is  15  kHz. 

The  two  DSPs  communicate  directly  through  the  dual-port  memory  bank  resident  on  the 
floating  point  DSP.  This  dual-port  memory,  which  is  divided  into  read-only  and  write- 
only  blocks  to  avoid  simultaneous  writes  to  the  same  location  by  the  two  processors, 
allows  the  two  processors  to  share  all  required  data  with  little  overhead.  The  floating 
point  DSP  provides  the  rotor  angle  estimate  and  current  commands  to  the  motor  control 
DSP,  while  the  motor  control  DSP  supplies  the  measured  currents  and  voltages  to  the 
floating  point  DSP. 


33  Software  Development 

In  the  DSP  implementation,  the  digital  controller  is  constructed  fiom  the  continuous-time 
observer  equations  (6),  (7),  (10)  and  (11)  using  the  forward  Euler  integration  method. 
The  development  of  the  software  was  carried  out  usii^  the  Analog  Devices  programming 
tools  for  the  Share  DSP,  and  the  naotion  control  debugger  for  the  fixed  point  motor 
control  DSP.  The  floating  and  fixed  point  devices  were  programmed  in  C  and  ^sembly 
language,  respectively.  The  PMSM  angular  velocity  and  angle  estimation  algorithms  are 
scheduled  to  execute  periodically  in  the  sequence  indicated  by  Figure  20,  which  shows 
program  flow  and  execution  times. 
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Current  Control  DSP  PMSM 


Figure  18.  Block  diagram  of  prototype  sensorless  drive  hardware. 


Figure  19.  Dual-DSP  processing  hardware-  top  view  (left),  bottom  view  (right). 

For  each  sampling  interval  of  the  floating  point  DSP  (100  ps),  the  observer  requires  input 
voltages  and  output  currents.  These  quantities  are  simultaneously  sampled  at  the  center 
point  of  each  PWM  cycle  [1]  by  the  motor  control  DSP  by  the  analog-to-digital  converter 
of  the  motor  control  DSP.  The  acquired  sensor  data  is  then  written  to  the  dual  port 
memory,  making  it  available  to  the  observer.  The  floating  point  DSP  operates  on  the 
acquired  data  to  update  the  estimated  state  variables  and  rotor  angle.  The  estimated  angle 
is  transmitted  back  to  the  motor  control  DSP  along  with  the  most  recent  velocity  estimate 
and  torque  command.  These  quantities  are  required  by  the  motor  control  DSP  to  execute 
the  torque  control  cycle.  The  floating  point  DSP  is  programmed  to  estimate  the  angular 
velocity  every  nth  cycle  time,  where  n=20  for  the  demonstrator  drive  system.  The 
observer  gam  matrix  G  is  updated  with  each  velocity  estimation  cycle. 
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Fixed  Point  DSP  Floating  Point  DSP 

A/D  Conversion  Complete  Timer  Interrupt  - 

Interrupt  (Center  of  PWM  Cycle)  Dual  Port  100  ^ 


Figure  20.  Program  flow  diagram  for  sensorless  drive  implementation. 


3.3.1  Sensors 

The  use  of  the  PWM  inverter  is  nearly  universal  in  PMSM  drives,  making  accurate 
voltage  and  current  measurements  problematic  due  to  the  switching  nature  of  the  applied 
voltage.  Due  to  the  electrically  noisy  environment  surrounding  a  typical  PWM  inverter, 
the  acquisition  of  high  quality  measurements  required  by  the  rotor  position  observer  is 
not  a  trivial  undertaking.  The  following  sections  describe  the  sensing  elements  and 
technique  used  to  acquire  the  PMSM  phase  vohages  and  currents. 


3  J.1.1  Voltage  Measurements 

The  approach  to  acquirii^  the  PMSM  phase  voltage  in  the  demonstrator  sensorless  drive 
is  the  direct  measurement  approach  [1].  Implementation  of  the  direct  measurement 
technique  requires  both  electrical  isolation  and  low-pass  filtering  in  order  to  isolate  power 
and  signal  grounds,  and  to  remove  high  fi'equency  switching  components  of  the  phase 
voltage,  re^wctively.  A  block  diagram  of  the  configuration  used  to  directly  measure  the 
phase  voltage  is  shown  in  Figure  21. 
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Figure  21 .  Block  diagram  of  direct  phase  voltage  measurement. 


In  Figure  21,  a  resistive  divider  is  used  to  scale  the  motor  terminal  voltage  and  to  form  a 
false  neutral  node,  which  under  normal  conditions  is  at  the  same  potential  as  the  buried 
PMSM  neutral.  The  voltage  transducer  (LEM  Model  CV  3-1000)  provides  an  isolated 
signal,  which  is  low-pass  filtered  before  being  sampled  by  the  motor  control  DSP  analog- 
to-digital  converter.  The  voltages  are  sampled  synchronously  to  the  PWM  fi-equency  at 
the  PWM  center  point,  as  described  in  [1].  This  sampling  strategy  reduces  the  effect  of 
PWM  ripple  on  the  sampled  data.  The  voltage  transducer  has  a  bandwidth  of  500  kHz  so 
that  it  does  not  distort  the  unfiltered  voltage  signal. 

The  low-pass  filter  is  a  4*  order  Bessel  filter  with  a  cutoff  fi-equency  of  5  kHz 
(Frequency  Devices  model  DP74L4L5kHz).  Despite  the  use  of  the  line  filter  between  the 
inverter  and  PMSM,  a  filter  is  still  required  to  further  remove  the  high  fi-equency 
switching  signal  fi-om  the  voltage  measurements.  The  drawback,  however,  is  that  the  use 
of  the  4*  order  filter  introduces  in  a  time  delay  of  approximately  68  pseconds  in  the 
voltage  signal.  This  delay  contributes  a  predictable  angle  estimation  error  at  higher 
fi-equencies.  This  phenomenon  will  be  shown  and  discussed  in  the  experimental  results 
of  Chapter  4. 

It  should  be  noted  that  the  indirect  voltage  measurement  technique  discussed  in  [1], 
which  uses  a  combination  of  the  DC  link  voltage  and  the  commanded  PWM  duty  cycles 
is  problematic  in  large  drive  systems  where  an  output  filter  is  used,  since  the  output  filter 
parameters  must  be  accurately  known  to  determine  the  voltages  present  at  the  motor 
terminal. 


3.3.1.2  Current  Measurements 

In  any  servo  application,  current  feedback  is  required  to  close  the  torque  control  loop. 
That  is,  current  measurements  are  required  in  either  sensor-based  or  sensorless  PMSM 
drives.  Like  the  voltage  signals,  the  ciurent  waveform  can  contain  higher  fi-equency 
components  due  to  the  PWM  driven  inverter.  Sampling  the  current  synchronously  with 
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the  PWM  waveform  relieves  this  situation  by  reducing  the  effect  of  PWM  ripple  on  the 
sampled  data.  In  doing  so  each  of  the  sideband  harmonics  are  symmetric  about  the  PWM 
switching  frequency  harmonic.  By  sampling  at  the  PWM  switching  frequency,  all  of  the 
sideband  harmonics  sum  to  DC.  Furthermore,  by  sampling  at  the  PWM  center  point,  the 
sideband  harmonics  sum  to  zero. 

The  demonstrator  sensorless  PMSM  drive  uses  200  ampere  current  transducers  (LEM 
Model  LA  205-S).  These  transducers  have  a  frequency  bandwidth  of  100  kHz  and  are 
very  linear  through  the  operating  range  of  the  device.  Only  two  of  these  sensors  are 
required  in  the  demonstrator  drive,  since  the  PMSM  is  connected  in  the  ‘  Y’  configuration 
and  the  third  phase  current  is  derived  from  the  two  measured  values.  Figure  22  shows  the 
current  sensors  as  installed  in  the  sensorless  drive  system  cabinet. 


Current  Sensors 


Figure  22.  Current  sensors  and  main  breaker. 


3.3.2  Power  Output  Stage 

The  power  output  stage  of  the  sensorless  electric  drive  system  serves  to  apply 
commanded  voltages  to  the  PMSM  terminals.  In  PMSM  drive  applications,  the  most 
popular  power  stage  consists  of  a  voltage  source  PWM-driven  inverter,  which  is  a  DC  to 
AC  converter  built  with  IGBTs  and  configured  to  deliver  bipolar  current  waveforms  to 
the  motor. 

In  the  demonstrator  sensorless  drive,  a  commercially  available  power  stage  is  used 
(Model  IPM1206,  Trace  Technologies,  Livermore,  CA).  This  inverter,  which  is  shown  in 
Figure  23,  is  capable  of  operating  at  DC  bus  voltages  of  up  to  750  volts,  and  (with  proper 
cooling)  is  capable  of  delivering  150  kW  at  PWM  switching  frequencies  of  up  to  15  kHz. 
The  DC  link  capacitor  bank  is  included  in  this  commercially  available  power  stage. 
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The  motor  control  DSP  provides  the  PWM  switching  signals  to  each  of  the  six  IGBT  gate 
drives  of  the  power  output  stage  through  a  fiber-optic  interface  (Agilent  Technologies 
Versalink  modules).  In  return,  the  power  stage  provides  fault  enunciation  for  each  IGBT 
through  the  same  fiber-optic  interface.  The  fiber-optic  interface  is  advantageous  due  to 
its  immunity  to  electrical  noise  and  inherent  electrical  isolation  of  power  electronics  and 


Figure  23.  Inverter  (left)  and  AC/DC  converter  (right)  used  in  sensorless  electric  drive. 
3.3.3  Output  Filter 

The  fast  voltage  risetime  delivered  by  modem  inverter  assemblies  (<1000  volts/ps)  has 
many  implications,  including  an  increase  m  common-mode  currents,  and  damaging 
overvoltage  conditions  at  the  motor  terminals  when  the  distance  between  the  motor  and 
the  inverter  is  long.  The  overvoltage  can  cause  early  breakdown  of  the  winding 
insulation  in  the  motor.  Common  mode  currents  may  result  in  increased  EMI  and 
increased  bearing  wear  due  to  the  introduction  of  shaft  currents. 

To  mitigate  the  undesirable  effects  of  the  steep  voltage  waveforms  present  in  high  power 
drive  assemblies,  a  filter  is  often  employed  to  limit  the  risetime  of  the  applied  voltage. 
The  output  filter  used  in  the  sensorless  electric  drive  is  comprised  of  a  single  stage  LC 
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filter,  as  shown  in  Figure  24.  The  LC  filter  as  shown  forms  a  low  pass  filter  with  a 
characteristic  2"*'  order  40  db/decade  attenuation  from  the  cutoff  frequency  of 


1 

2;rVZc  ■ 


(78) 


The  LC  filter  used  (MTE  Corporation  model  RL10013C)  is  a  75  HP  unit  that  combines 
the  characteristics  of  line  reactors,  which  are  known  to  limit  common  mode  currents,  with 
the  capacitor  network  to  generate  a  near-sinusoidal  output  voltage  waveform,  as  shown  in 
Figure  25.  The  use  of  such  filters  is  seen  to  significantly  reduce  the  dv/dt  of  the  output 
voltage  and  thus  helps  to  extend  motor  life  by  reducing  the  motor  peak  voltage. 


Figure  25  Line  voltage  without  LC  filter  (left)  and  with  LC  filter  (right). 

3.4  PMSM  Motor 

The  experimental  motor  used  in  the  sensorless  electric  drive  demonstrator  is  a  70  HP 
PMSM  manufactured  by  Electrodynamics  Corp  division  of  Electric  Boat,  Avanel  NJ. 
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The  machine  was  designed  for  dual-use,  and  therefore  contains  separate  phase  winding 
connections,  as  depicted  in  Figure  26.  In  the  low  speed  configuration,  the  individual 
windings  in  each  of  the  three  phases  are  connected  in  series,  while  for  high  speed 
applications,  they  are  connected  in  parallel.  The  parameters  of  the  machine  are  shown  in 
Tables  I  and  II  for  the  low  and  high  speed  configurations,  respectively. 


Figure  27.  70  HP  PMSM  used  in  sensorless  electric  drive  studies. 
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Table  I.  Electrodynamics  Corp.  motor  parameters  (low-speed  winding  configuration). 


Parameter 

Value 

Inductance  Lav 

1.7  mH 

Resistance  R 

.28  Q 

PM  Flux  Linkage 

1.376  v-s 

Rated  Power 

70  HP 

Poles 

8 

Table  II.  Electrodynamics  Corp.  motor  parameters  (high-speed  winding  configuration). 


Parameter 

Value 

Inductance  Lav 

.425  mH 

Resistance  R 

.07  Q 

PM  Flux  Linkage 

0.688  v-s 

Rated  Power 

70  HP 

Poles 

8 

4  SENSORLESS  PMSM  EXPERIMENTAL  RESULTS 

In  this  chapter,  the  rotor  position  estimation  strategy  is  operated  in  a  closed-loop  drive 
system,  where  the  drive  is  implemented  as  discussed  in  Chapter  3.  To  evaluate  the 
performance  of  the  sensorless  drive,  a  set  of  experiments  were  performed  where  the 
computing  hardware  executes  the  complete  PMSM  control  algorithm  in  real  time. 
Experimental  results  depicting  the  rotor  position  estimation  accuracy  across  the  operating 
range  of  the  PMSM  are  presented  in  this  chapter  for  both  transient,  steady-state,  and 
varying  load  torque  conditions. 


4.1  Experimental  System 

The  test  setup  used  in  the  experiments  is  shown  in  Figure  1 1  (drive  electronics)  and 
Figure  28  (motor/load).  The  mechanical  arrangement  consists  of  a  PMSM  whose  shaft 
is  coupled  to  a  pump,  which  is  used  to  provide  load  torque.  A  quadrature  shaft  encoder 
(14,400  counts  per  mechanical  revolution)  was  attached  directly  to  the  motor  shaft.  The 
encoder  was  not  employed  by  the  drive  system,  but  was  used  to  measure  the  accuracy  of 
the  estimated  rotor  angle.  The  drive  electronics  shown  in  Figure  1 1  and  the  motor  shown 
in  Figure  28  are  physically  separated  by  approximately  75  feet  of  cable. 
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Figure  28.  Motor  (foreground)  and  load  (background)  used  in  sensorless  drive  studies. 
4.2  Experimental  Results 

In  this  section,  experiments  on  the  sensorless  drive  are  carried  out.  All  of  the 
experiments  are  executed  using  the  developed  sensorless  drive  hardware,  which  controls 
the  PMSM  in  real-time  using  the  estimated  rotor  angle  and  speed.  The  experiments  are 
intended  to  demonstrate  the  effectiveness  of  the  proposed  sensorless  technique  over  a 
wide  variety  of  operating  conditions. 


4.2.1  Steady  State  Rotor  Position  Estimation  Accuracy 

To  determine  the  angle  estimation  accuracy,  a  set  of  experiments  was  performed  on  the 
drive  system  across  the  speed  range  of  the  PMSM.  The  results  are  shown  in  Figure  29 
and  Figure  30,  where  the  rotor  position  estimation  error  is  shown  for  the  low  speed  and 
high  speed  winding  configurations,  respectively.  For  these  tests,  the  PMSM  was 
essentially  unloaded.  It  will  be  demonstrated  in  a  later  section  that  the  angle  estimation 
error  is  largely  independent  of  load  torque. 
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Angular  Velocity  (RPM) 


Figure  29.  Ai^le  estimation  accuracy  versus  speed  for  low  speed  winding  configuratioa 

RPM  vs.  Angle  Est.  Error  •  High  Speed  Windings 


Angular  Velocity  (RPM) 


Figure  30.  Angle  estimation  accuracy  versus  speed  for  high  speed  winding  configuration. 


Although  the  performance  of  the  rotor  position  estimator  is  excellent  in  Figure  29  and 
Figure  30,  it  is  clear  that  the  estimation  error  increases  with  rotor  speed  (frequency).  This 
degradation  in  rotor  position  estimation  performance  as  speed  increases  is  due  to  the 
delay  introduced  by  the  low-pass  filtering  of  the  phase  voltages  and  the  sampling  delays 
introduced  by  the  discrete-time  implementation  of  the  algorithm-  The  filter  and  sampling 
delays  are  67  ps  and  100  ps,  respectively.  Thus,  the  net  effect  is  a  linear  phase  shift  of 
the  measured  phase  voltages  that  is  linearly  dependent  on  the  operating  frequency.  The 
phase  shift  produces  a  corresponding  frequency-dependent  angle  estimation  error  given 

by 


phase  shift  (r)  =  167/is  *  Shaft  Speed  (in  r/ 5)  .  (79) 

For  convenience,  this  expected  angle  estimation  error  due  to  the  filtering  and  sampling 
delays  given  in  (79)  is  included  in  Figure  29  and  Figure  30.  Clearly,  the  angle  e^unation 
error  follows  the  expected  error  due  to  the  delays.  It  should  be  noted  that  this  speed- 
dependent  estimation  error  could  be  easily  removed  by  compensating  for  the 
filter/sampling  delays,  since  the  error  versus  speed  characteristic  is  well  known  from 

(79). 


4.2.1.1  Low  Speed  Operation 

In  the  sensorless  drive,  rotor  position  estimation  at  and  around  zero  speed  present  is 
problematic.  In  the  demonstrator  system,  this  is  managed  by  setting  the  observer  gain 
matrbc  G  to  zero  for  angular  velocities  below  an  experimentally  determined  low-speed 
threshold,  which  essentially  disables  the  observer.  Under  these  operating  conditions,  the 
internal  PMSM  model  runs  without  correctional  feedback.  The  low-speed  threshold  is 
7.5  r/s,  or  approximately  19  RPM  for  the  8  pole  PMSM  used  in  the  studies. 

Figure  31  and  Figure  32  show  the  operation  of  the  sensorless  drive  at  extremely  low 
speeds.  Figure  31  demonstrates  operation  at  a  speed  slightly  above  the  low-speed 
threshold  of  19  RPM.  In  this  experiment,  the  motor  runs  quite  reliably.  However, 
occasional  polarity  reversals  in  the  estimated  velocity  occur  at  such  low  speeds  that  tend 
to  generate  correspondii^  deviations  in  the  estimate  rotor  angle.  In  situations  such  as 
Figure  31,  where  polarity  reversals  are  rather  infrequent  and  of  short  duration,  the 
estimated  rotor  angle  error  diverges  during  the  estimated  velocity  polarity  reversals,  but 
is  quickly  adjusted  when  the  angular  velocity  estimate  is  corrected. 

In  Figure  32,  the  system  is  operating  at  50  RPM,  which  is  more  than  twice  the  low-speed 
threshold.  At  this  speed,  the  estimated  velocity  polarity  reversals  do  not  occur,  and  the 
angle  estimation  error  is  much  more  stable.  From  a  comparison  of  Figiure  3 1  and  Figure 
32,  it  is  apparent  that  the  estimation  of  angular  velocity  is  more  problematic  at  lower 
sp^s.  This  is  due  to  the  difficulty  in  determining  the  direction  of  rotation  of  the  voltage 
vector  described  in  (15)  and  (16)  when  noise  is  present  on  the  low-amplitude  terminal 
voltage  signals. 
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Figure  31,  Estimation  accuracy  at  30  RPM,  10  ft-lb  -  low  speed  winding  configuratioa 


Figure  32.  Estimation  accuracy  at  50  RPM,  15  ft-lb  -  low  speed  winding  configuration. 
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For  the  sensorless  drive  system  under  consideration,  the  low-speed  threshold  was 
experimentally  determined  by  viewing  the  estimated  angular  velocity  at  several  low 
speeds.  The  speed  at  which  the  polarity  reversals  become  quite  infrequent  represents  a 
reasonable  selection  for  the  low-speed  threshold.  It  has  been  the  experience  of  the  author 
that  this  low-speed  threshold  is  typically  2%  to  4%  of  the  rated  motor  speed. 


4.2.1.2  Steady  State  Rotor  Position  Estimation  Robustness 

In  this  section,  the  robustness  of  the  developed  sensorless  PMSM  drive  is  examined  with 
respect  to  parameter  uncertainty  and  external  disturbances.  In  practical  applications, 
machine  parameters  vary  due  to  temperature  effects.  Therefore,  an  effective  rotor 
position  estimator  must  be  robust  to  any  uncertainty  in  the  machine  model. 

The  effect  of  load  torque  disturbances  is  shown  in  Figure  33  and  Figure  34  for  the  low 
speed  and  high  speed  winding  configuration,  respectively.  Recall  that  the  model  given  in 
(6)  and  (7)  does  not  contain  mechanical  variables  such  as  torque,  friction,  and  inertia  in 
the  state  equations.  Thus  variations  in  the  load  torque  are  not  expected  to  significantly 
affect  the  accuracy  of  the  rotor  position  estimator.  However,  it  should  be  noted  that  an 
increase  in  load  torque  requires  a  corresponding  rise  in  the  stator  current.  From  the  ai^le 
estimation  error  analysis  of  Chapter  2,  it  can  be  seen  that  a  torque  increase  can  cause  a 
change  in  angle  estimation  accuracy  if  there  is  uncertainty  in  the  motor  inductance  or 
resistance  parameters,  or  if  there  exists  errors  in  the  current  measurement.  The  predicted 
and  experimentally  certified  immunity  of  rotor  position  estimation  accuracy  to  load 
torque  is  verified  in  both  Figure  33  and  Figure  34.  This  represents  a  tremendous 
advantage  over  many  other  sensorless  drive  strategies  that  require  estimation  or  direct 
measurement  of  the  load  torque. 

The  effect  of  stator  inductance  uncertainty  on  position  estimation  accuracy  has  been 
examined  in  Chapter  2.  In  the  experiment  depicted  in  Figure  35,  the  stator  inductance 
used  by  the  observer  PMSM  model  is  intentionally  detuned  by  50%  from  its  actual  value. 
Comparison  of  Figure  33  and  Figure  35  show  how  the  performance  of  the  rotor  position 
estimation  accuracy  is  adversely  affected  by  the  changing  load  torque  in  the  case  of  stator 
inductance  uncertainty.  Clearly,  the  inaccurate  inductance  used  by  the  model  results  in 
degraded  performance  with  increasing  load  (current)  as  predicted  in  Chapter  2. 
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Load  Torqua  (Ft-Lb) 


Steady  State  Angle  Estimation  -  Varying  Load 


Figure  35.  Estimation  accuracy  versus  load  at  Constant  Speed  -  low  speed  windings 
with  stator  inductance  detuned  by  50%. 


In  the  PMSM,  the  machine  parameters  and  R  may  vary  inversely  and  proportionally 
with  the  temperature,  respectively.  The  permanent  magnet  flux  linkage  parameter  is 
not  used  in  the  machine  model  so  that  it  is  not  expected  to  directly  generate  angle 
estimation  errors.  However,  the  parameter  is  used  to  estimate  the  angular  velocity  so 

that  uncertainty  in  A„  manifests  itself  as  a  velocity  estimation  error.  Note,  however,  that 
even  with  the  adaptive  angular  velocity  estimation  technique  presented  in  Chapter  2, 
reasonable  uncertainty  in  A„  will  not  produce  angular  velocity  errors  in  the  steady  state. 

In  Chapter  2,  the  effects  of  stator  resistance  uncertainty  were  examined.  This  study 
revealed  that  under  normal  operating  conditions,  stator  resistance  uncertainty  has  little 
effect  on  angle  estimation  accuracy.  To  verify  this  prediction,  angle  estimation  error 
measurements  were  performed  with  the  internal  model  parameters  intentionally  detuned 
by  (an  extreme)  300  percent.  The  results  of  this  experiment  are  shown  in  Figure  36, 
where  the  top  plots  show  performance  with  the  actual  resistance  value.  The  bottom  plots 
show  system  performance  when  the  PMSM  model  uses  a  resistance  value  that  is 
increas^  by  300%  from  its  known  value.  In  both  cases,  the  motor  was  operating  at  200 
RPM  with  a  load  of  60  fl-lb.  From  a  comparison  of  the  results  in  Figtire  36,  the 
prediction  set  forth  in  the  estimation  error  analysis  of  Chapter  2  is  confirmed.  It  should 
be  noted  that  the  resistance  parameter  that  is  intentionally  being  detuned  is  also  used  to 
determine  the  values  in  the  observer  gain  matrix  G  given  in  (23),  and  to  determine  the 
estimated  angular  velocity  that  is  used  to  construct  the  time- varying  system  matrk  (Ao,). 
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Figure  36.  Estimation  accuracy  with  actual  (top),  and  detuned  resistance  values. 


One  concern  in  the  implementation  of  a  sensorless  electric  drive  is  related  to  the  use  of 
the  line  filter  between  the  inverter  and  PMSM,  as  shown  in  Figure  24.  When  such  a  filter 
is  utilized,  it  is  important  to  sense  the  voltage  and  current  at  the  filter  output  as  shown  in 
Figure  10,  rather  than  at  the  inverter  output.  For  this  reason,  the  line  filter  presents 
difficulties  in  using  the  indirect  voltage  sensing  method,  as  described  in  [18]. 

To  examine  the  effects  of  the  line  filter  on  rotor  position  estimation,  a  set  of  experiments 
were  carried  out.  The  drive  was  operated  at  various  RPM  and  load  conditions  both  with 
and  without  the  line  filter  in  place,  and  the  results  were  conqjared.  An  example  of  the 
test  results  is  presented  in  Figure  37  and  Figure  38,  where  the  PMSM  is  operating  at  200 
RPM  and  100  ft- lb  load  torque.  Figure  37  shows  the  performance  with  the  filter  in  place, 
while  performance  without  the  filter  is  shown  in  Figure  38.  Comparison  of  these 
experimental  results  demonstrates  little  if  any  change  in  performance.  This  result  was 
confirmed  by  the  test  results  at  other  operating  conditions  across  the  range  of  the  PMSM 
capabilities. 

A  complete  plot  of  the  measured  quantities,  estimated  states,  and  estimated  rotor  angle  is 
included  in  Figure  39  and  Figme  40.  Figure  39  is  shown  for  the  low-speed  winding 
configuration  operating  at  500  RPM  and  400  ft-lb.  load  torque.  Figure  40  shows 
operation  at  900  RPM  and  400  ft-lb  (70  HP).  In  each  of  these  figures,  note  that  the 
measured  quantities  (voltage  and  current)  are  shown  in  the  stationary  three-phase 
reference  frame.  All  estimated  values  (flux  linkages  and  currents)  are  shown  in  the  two- 
phase  stationary  reference  flume. 
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4.2.2  Transient  Rotor  Position  Estimation  Accuracy 

During  startiq)  or  speed  reversal  transients,  the  rotor  velocity  passes  through  the 
problematic  zero  speed  range.  Therefore,  an  examination  of  the  transient  performance  of 
the  sensorless  drive  is  essential.  In  this  section,  the  system  performance  is  examined  for 
several  transient  conditions,  including  startup  and  speed  reversal. 


4.2.2. 1  Startup  of  Sensorless  PMSM 

It  is  well  known  that  sensorless  operation  is  problematic  at  standstill.  At  zero  speed, 
general  purpose  principles  for  rotor  position  estimation  usually  fail  due  to  the  lack  of  an 
induced  voltage  in  the  stator  windings.  For  the  proposed  rotor  position  estimator,  the 
observer  is  ill-conditioned  at  zero  speed  because  the  last  two  rows  of  the  system  matrix  A 
go  to  zero,  rendering  the  system  unobservable.  Therefore,  the  position  estimator  is  not 
usable  at  zero  speed,  and  unless  the  initial  rotor  position  is  known,  it  is  not  possible  to 
start  the  machine  at  full  torque  or  to  guarantee  a  dither-free  startup. 

For  the  sensorless  drive  demonstrator  system,  a  specialized  startup  strategy  is 
implemented  that  yields  to  the  rotor  position  observer  when  a  low-speed  threshold  is 
exceeded.  More  specifically,  the  demonstrator  sensorless  drive  adopts  the  so-called 
forced-alignment  technique  [24].  The  forced  alignment  technique  ^lies  a  dc  current  to 
the  stator  before  startup.  The  dc  current  acts  to  align  the  permanent  magnet  field  with  the 
magnetic  field  generated  by  the  stator  excitation,  and  forces  the  rotor  to  a  known  initial 
position.  Since  the  initial  position  is  known,  torque  can  be  applied  to  the  PMSM  at 
startup  and  the  rotor  position  estimator  operated  in  open  loop  configuration  until  the  low- 
speed  threshold  is  surpassed. 

The  initial  torque  applied  to  the  PMSM  at  standstill  after  forced  alignment  must  be 
sufficient  to  produce  a  rotational  velocity  that  exceeds  the  low-speed  threshold  before  the 
angle  estimation  error  becomes  excessive.  In  practice,  this  objective  has  been  achieved 
to  produce  very  reliable  PMSM  startup  despite  the  lack  of  a  position  sensor.  To 
demonstrate  the  startup  procedure  (after  forced  alignment)  several  illustrations  are 
provided.  In  the  startup  depicted  in  Figure  41  where  a  startup  velocity  command  of  200 
RPM  is  issued,  the  low-speed  threshold  of  7.5  r/s  is  exceeded  at  0.96  seconds.  At  that 
instant,  the  estimated  rotor  angle  converges  toward  zero  fix)m  an  initial  error  of 
approximately  20  electrical  degrees.  Until  the  low-speed  threshold  is  achieved,  the 
system  operates  without  feedback,  with  the  estimated  rotor  velocity  driving  the  harmonic 
oscillator  that  produces  the  estimated  permanent  magnet  flux  linkages  in  the  two-phase 
stationary  reference  frame. 

In  Figure  42,  each  of  the  measured  quantities  (current  and  voltage)  as  well  as  the 
estimated  state  variables  are  shown  for  a  startup  sequence.  For  this  experiment,  a  speed 
command  of  300  RPM  is  issued,  with  the  viscous  load  set  to  provide  1 80  ft-lb  of  torque 
at  the  commanded  speed. 
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startup  -  80  Ft-Lb  at  200  RPM 


Figiire  41 .  Startup  from  zero  speed  with  load  of  80  fl-lb  at  200  RPM. 


4.2.2.2  Speed  Reversal 

During  a  speed  reversal,  the  rotor  velocity  passes  through  the  problematic  zero  speed 
range.  While  the  magnitude  of  the  estimated  rotor  velocity  is  below  the  low-speed 
threshold,  the  rotor  position  observer  operates  without  feedback.  Such  operation  has 
been  shown  through  experimentation  to  be  acceptable  for  short  time  intervals,  such  as 
during  a  speed  reversal  In  Figure  43,  the  rotor  position  observer  is  examined  during  a 
speed  reversal  In  the  figure,  the  estimated  velocity  magnitude  drops  below  the  threshold 
at  1.75  sec.  And  the  observer  operates  without  feedback  until  1.9  seconds,  where 
feedback  is  restored.  From  the  estimated  states  in  the  figure  (flux  linkage  and 
trigonometric  functions  of  the  rotor  angle)  it  is  seen  that  the  open- loop  model  produces 
excellent  results,  or  ride-through  capability,  for  short  term  operation  below  the  low-speed 
threshold. 
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43  Acoustic  Performance 


To  quantify  the  effects  of  the  sensorless  drive  on  acoustic  performance,  a  set  of  self-noise 
measurements  are  to  be  performed  on  the  drive,  and  compared  with  sensored  operation  of 
the  same  motor.  As  indicated  in  the  original  proposal,  these  tests  will  be  provided  to 
ONR  for  at  least  the  first  two  water  tunnel  tests  that  use  the  developed  sensorless  drive 
system.  At  the  time  of  this  writing,  there  have  not  been  any  water  txmnel  tests  scheduled 
to  utilize  the  PMSM  and  thus  the  data  is  not  yet  available.  An  addendum  to  this 
document  will  be  generated  at  the  time  this  acoustic  measurement  data  becomes  is 
obtained.  As  indicated  in  the  proposal,  the  cost  for  these  tests  will  be  absorbed  by  ARL. 


5  CONCLUSIONS 

This  work  has  presented  the  design  and  implementation  of  a  sensorless  electric  drive  for  a 
high  power  permanent  magnet  synchronous  motor  (PMSM).  The  developed  sensorless 
drive  is  implemented  in  a  real-time  nootor  drive  system,  and  is  shown  through  extensive 
experimentation  to  be  robust  to  model  uncertainty  and  disturbances,  while  providing  very 
accurate  rotor  position  estimation.  The  angle  estimation  accuracy  demonstrated  by  the 
proposed  technique  exceeds  the  published  accuracy  data  fiom  other  laboratory  sensorless 
drive  implementations.  Furthermore,  the  proposed  technique  has  been  shown  to  be 
itniniine  to  mechanical  parameters  and  disturbances,  which  are  difficuh  to  determine  and 
are  often  unknown. 

The  primary  remaining  issue  associated  with  the  use  of  the  sensorless  drive  in  a 
propulsion  system  is  the  implementation  of  a  high-performance  velocity  controUer.  Due 
to  the  high  accuracy  demonstrated  by  the  rotor  position  estimator,  current  control 
performance  approaches  that  of  a  fully  sensored  drive,  with  clear  advantages  of  sensor- 
based  control  only  evident  aroimd  zero  speed.  The  lack  of  a  velocity  sensor  in  the 
sensorless  system,  however,  presents  a  serious  challenge  in  achieving  accurate  velocity 
tracking.  As  discussed  in  tte  text,  the  adaptive  velocity  estimation  technique  provides 
accurate  velocity  feedback  in  the  steady  state,  but  low-pass  fikering  with  a  very  slow 
time  constant  is  required  to  remove  ripple  from  the  estimated  angular  velocity.  If  the 
ripple  is  not  remov^  the  torqiw  control  loop  suffers,  and  smooth  torque  is  not  generated. 
Excessive  filtering  of  the  estimated  velocity,  however,  degrades  the  transient 
performance  of  the  system. 

Issues  related  to  sensorless  PMSM  startup  have  been  addressed  in  this  document  through 
the  adoption  of  a  forced  alignment  tech^que.  This  technique  has  been  shown  through 
experimentation  to  provide  reliable  PMSM  startup. 

In  summary,  the  results  of  this  study  suggest  that  the  sensorless  electric  drive  is  a  viable 
alternative  for  electric  drive  applications  such  as  rim-driven  propulsion  topologies.  The 
sensorless  drive  i^^)lementation  and  experimental  results  discussed  in  this  study  will 
serve  as  valuable  information  for  any  future  sensorless  drive  design  endeavors. 
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